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Abstract. Discoveries of very high energy (VHE) photons from distant blazars suggest 
that, after correction by extragalactic background light (EBL) absorption, there is a flatness 
or even a turn-up in their spectra at the highest energies that cannot be easily explained 
by the standard framework. Here, it is shown that a possible solution to this problem is 
achieved by assuming the existence of axion-like particles (ALPs) with masses ~ 1 neV. The 
ALP scenario is tested making use of observations of the highest redshift blazars known in the 
VHE energy regime, namely 3C 279, 3C 66A, PKS 1222+216 and PG 1553+113. In all cases, 
better fits to the observed spectra are found when including ALPs rather than considering 
EBL only. Interestingly, quite similar critical energies for photon/ ALP conversions are also 
derived, independently of the source considered. 
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1 Introduction 

Blazars are active galactic nuclei (AGN) characterized by having their jets oriented at small 
angles towards our line of sight. Their emission covers the entire electromagnetic spectrum 
from radio up to the most energetic 7-rays. The synchrotron self Compton (SSC) model 
provides a successful description for the double-peaked broad spectral energy distribution of 
low redshift blazars up to several TeV e.g., [1, 2]. 

Very high energy photons (VHE, 30 GeV-30 TeV) coming from cosmological sources, 
such as blazars, are attenuated by pair-production with ultraviolet (UV), optical, and infrared 
(IR) photons in the extragalactic background light (EBL). The EBL is the accumulated 
radiation due to star formation plus a contribution from AGN emission in galaxies over the 
entire history of the Universe. As summarized in refs. [3, 4], the EBL coming from galaxies 
seems already well constrained from UV up to mid-IR wavelengths at the minimum intensity 
level allowed by galaxy counts. This implies the highest transparency to 7-ray photons from 
standard physics. Independent EBL data from galaxy counts [5-7] and different modelings 
are already in good agreement [3, 8-10], even though direct detection remains elusive due 
to huge contamination from the zodiacal light [11]. However, a very recent analysis of a 
data set, which the authors claim to be zodiacal-light free, seems to support the EBL at an 
intensity level close to galaxy counts as already suggested by other methods [12]. On the 
other hand, discoveries of VHE photons from high redshift sources by imaging atmospheric 
Cherenkov telescopes (IACTs) such as MAGIC, VERITAS, and HESS also suggest an EBL 
intensity level in agreement with galaxy counts and EBL models, e.g., [13-16]. Indeed, 
at present, the uncertainties in the recovered unattenuated VHE spectra (hereafter, EBL- 
corrected spectra) of blazars observed with the present generation of IACTs are in general 
dominated by statistical and systematic effects rather than uncertainties in the EBL modeling 



When referring to differential fluxes, EBL-corrected spectra of some blazars suggest a 
flatness or even an upturn at the highest observed energies (hereafter, the pile-up problem; see 
e.g., figure 3 in [4]), which cannot be easily described by the standard framework described 
above. The pile-up problem was already noticed some years ago (see e.g., refs. [17, 18]) 
when a higher intensity EBL was expected. At the time, the most possible solution seemed 
to be to assume lower EBL intensities (as done e.g., in ref. [13]). However, as already 



[3, 4]. 
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discussed, today independent EBL results seem to converge towards a Universe with the 
highest possible transparency to 7-ray photons. Therefore, the pile-up problem, if present, 
can in general no longer be attributed to EBL issues. Modifications to the so-called standard 
one- zone SSC model have also been proposed in order to solve very hard VHE spectra such as 
strong relativistic shocks [19], selection bias effects [20], internal absorption [21], production 
in extended jets [22], but there is still little work focused in solving the pile-up problem 
e.g., [17]. Indeed, there is not any general agreement on any definitive solution yet. The 
postulated existence of axion-like particles (ALPs) arises as an alternative solution to the 
problem. In fact, the main goal of this paper is to show explicitly that the theoretical 
framework developed in ref. [23] when applied to measured VHE spectra of distant blazars 
is a possible solution that may deserve further attention. 

2 Theoretical framework 

Axions are predicted by the Peccei-Quinn mechanism, which is currently the most compelling 
explanation to solve the CP problem in quantum chromo dynamics [25]. They might consti- 
tute a portion or the totality of the non-baryonic cold dark matter content in the Universe 
[26]. One may also consider the existence of ALPs, which are particles with the same prop- 
erties as axions, but with mass and coupling constant not related to each other [27]. There is 
an additional property of ALPs that could have important implications for their detectabil- 
ity as well as for 7-ray astronomy, i.e., they can convert into photons and vice- versa in the 
presence of an electric or magnetic field [28, 29]. This argument was first investigated in the 
optical band by ref. [30], where authors proposed the existence of axions to be the cause of 
the observed supernova la dimming (this proposal is ruled out at present; see e.g., ref. [31] 
and references therein). Later, the photon/ ALP conversion was applied for the same authors 
at higher energies in ref. [32]. The complementary view of photon/ ALP oscillations both in 
the source and in the Galactic magnetic fields was first put forward in ref. [33]. It has been 
also proposed that ALPs might significantly affect 7-ray AGN observations [34-37] . Ref. [23] 
was the first work in offering a complete formalism that includes under the same consistent 
framework the photon/ ALP mixing expected to occur in both the 7-ray sources and in the 
intergalactic medium (IGM). There, prospects for ALP detection were also discussed. One 
possible approach lies in detecting intensity boost in the higher energy bins of cosmological 
VHE observations by IACTs as a consequence of an efficient photon/ ALP conversion. Other 
approaches might be considered as well, such as detecting abrupt changes in intensity at a 
given critical energy, E cr i t , defined in GeV as E cr u = m^yMn/OABc, which should be the 
same for all the sources [34, 38]. In this equation, M%\ is the photon/ALP coupling constant 
divided by 10 11 GeV, B G the magnetic field in units of G, and m^ e y the effective ALP mass in 
^eV, defined as = \rr? a — a^J, where m a is the ALP mass and oj p i = 0.37 x lQ" A y /n e the 

plasma frequency in [ieV, n e being the electron density in units of cm -3 . These alternative 
approaches will be followed elsewhere. 

3 Application to distant known blazars 

Our methodology for studying any possible hint of a signature introduced by ALP in the 
observed VHE blazar spectra is based on several assumptions: (i) the intrinsic spectra emitted 
from the sources are well described by simple power laws, which is a good approximation 
for the relatively small energy ranges considered in the blazars studied here; (ii) M\\ is an 
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Figure 1. The effect of ALPs on the spectrum of 3C 279 obtained from 2006 MAGIC data [15]. Data 
points corrected for EBL only and for EBL plus ALP in Steps 1 and 2 are slightly shifted for clarity. 
See the text for a detailed explanation of every step. 



optimistic value, so that we can get maximum photon/ALP conversions, but still within 
present experimental limits, i.e., those imposed by the CAST experiment [39] 1 ; (iii) the value 
of the critical energy corresponding to the intergalactic mixing (hereafter, E cr n) lies in the 
energy range measurable by IACTs, or below; (iv) the EBL is well described by the model 
described in ref. [3]. 

The effect of ALP on 7-rays being emitted by blazars is modeled as follows: (i) the 
intrinsic emission from every source is obtained from the best SSC model fit to the multi- 
wavelength observations for each source; (ii) through the IGM we adopt two different values 
of the intergalactic magnetic field Bigm = 0.1 nG and Bigm = 1 nG given its large uncer- 
tainties at present e.g., [40]; (iii) once E cr n and Bigm are fixed, we then calculate the ALP 
mass m a using the relation between these quantities described in the above formula. In this 
study, we scan E cr n within the range 50-625 GeV in steps of 25 GeV, which corresponds to 
m a ~ 1 neV. It is important to note that any E cr u below 50 GeV will lead, in the energy 
range under consideration, to identical results than those obtained for 50 GeV. The reason 
is the strong dependence of the ALP effect on the EBL intensity. As the role of the EBL 
is negligible below ~ 50 GeV, photon/ALP mixings up to this energy will produce identical 
features in the spectra. Relevant source, IGM, and ALP parameters are listed in Table 1. 

We now carefully describe our methodology when applied to 3C 279, which is later 
repeated for other potentially interesting blazars, in particular 3C 66A, PKS 1222+216 and 
PG 1553+113. These four objects were selected according to their high redshift, where the 
effect of photon/ALP conversion in observed spectra is expected to be larger at the highest 
measured energies [23]. Indeed, this sample actually contains the four highest-redshift blazars 
observed so far and for which there is published spectra available. 

3.1 3C 279 

We show in figure 1 and figure 2 (upper-left panel) VHE data from 3C 279 [z = 0.536). 
These spectra were measured by the MAGIC collaboration at two different epochs [15, 41] 
when this flat spectrum radio quasar (FSRQ) was in different flaring states. To model the 

1 We note that significantly higher M values will probably lead to a negligible effect due to photon/ALP 
conversions in the sources. Yet, the effect coming from conversions in the intergalactic magnetic field is not 
negligible, but would be substantially lower than the one shown in this paper. As a consequence, ALPs might 
not probably work fine to solve the pile-up problem in many sources. 
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Figure 2. The effect of ALPs on high redshift VHE observations. Upper-left panel: 3C 279 obser- 
vation by MAGIC in the 2007 observational campaign [41]. Upper-right panel: 3C 66A observation 
by MAGIC [42]. Lower-left panel: PKS 1222+216 discovery by MAGIC [16]. Lower-right panel: 
PG 1553+113 extensive campaign from 2005-2009 by MAGIC [43]. As additional case, the 3C 279 
spectrum observed by MAGIC in 2006 was already shown in figure 1 (Step 4). 



photon/ALP mixing within the source, we use the SSC parameters given for both epochs in 
ref. [41] (see Table 1). Different steps are followed in order to obtain our final results. These 
steps are only explicitly shown in figure 1 for the 2006 data, while only the last step is plotted 
in figure 2 for the other spectra used in this study. Step 1) The observed VHE spectrum 
(blue crosses) is corrected by the effect from the EBL only (green circles), and by the EBL 
plus ALPs (red squares). The case with E CT n = 100 GeV and Bigm = 0.1 nG is plotted in 
figure 1 as an example. The pile-up problem is evident here for the EBL-corrected spectrum 
(green circles). Step 2) Both corrected (or de-absorbed) spectra are fitted by simple power 
laws. The spectral indexes and the Xr = X 2 l n °f the fits are shown in figure 1 (with n 
degrees of freedom). Step 3) The best-fit power laws obtained from Step 2 are then absorbed 
again by the EBL (dashed green line), and the EBL plus ALP (solid red line). Although 
small, the drop in flux at E cr n (in this particular case 100 GeV) is still noticeable. This step 
is taken in order to check how well the observed spectrum is recovered in the case where 
we account for only the EBL effect and when including EBL plus ALPs. The reduced \ 2 
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labeled in the figure corresponds to the corrected best-fit power law given in Step 2. Step 
4) Steps 1 to 3 are then repeated by varying E cr u in 25 GeV steps within the energy range 
in consideration and using the two different values of Bjgm considered here (see Table 1). 
The best fit among all the cases is plotted in figure 1 (right panel) as a solid red line. For 
this particular source, the best fit is obtained for E cr a = 100 GeV. Note that this is indeed a 
better fit to the data than when only considering the EBL attenuation (dashed green line). 
This is quantitatively shown in Table 2, where we provide, for all sources, the Xr °f the fits 
for five benchmark E cr u in the spectra. Finally, the shadow region in every figure represents 
the set of fits obtained when varying E cr it within the whole energy interval, in 25 GeV steps 
as mentioned above. Therefore, this shadow region might be taken as the uncertainty region 
that arises due to our lack of knowledge of the exact E cr u. Note that E cr u does not enter as 
an additional parameter in the fits; indeed, both kind of fits (i.e., with and without ALPs) 
do actually have the same degrees of freedom. The uncertainty region simply reflects a lack 
of knowledge on the underlying physics (real ALP parameters and Biqm), each particular 
fit in this region being completely independent of the others (meaning that E cr a does not 
act MS M linking parameter between them). 

3.2 3C 66A 

This source at z = 0.444 was observed in a flaring state by VERITAS [44] and more recently 
by MAGIC [42]. Both spectra are roughly compatible within uncertainties. Note that the 
redshift considered for this object is calculated using just one emission line and is thus not 
very secure, as discussed in ref. [45]. In this study, only the MAGIC spectrum is used due to 
its larger energy coverage, which better exposes the pile-up problem. We model the source 
using the parameters given in Table 1 , which were extracted from ref. [46] and were computed 
by the authors using an SSC model. Figure 2 shows the 3C 66A spectrum as measured by 
MAGIC (upper-right panel). The best fit is shown following the same methodology as in 
figure 1 for the case of 3C 279. We only show the results corresponding to Step 4. Again, 
considering the existence of ALP provides a better fit. 

3.3 PKS 1222+216 

This FSRQ located at z = 0.432 was discovered in the VHE regime in a flaring state by the 
MAGIC collaboration [16]. We model the source using the SSC parameters given in [47] for 
the so-called all blob emission model. In any case, we stress that the effect of the photon/ALP 
conversion within the source is estimated to be less than 10% of the total. In this particular 
case, the mixing within the source is negligible mainly due to the small value of the region 
where the magnetic field is confined (R, see Table 1) and therefore only intergalactic mixing 
is present. The results obtained for this blazar are shown in figure 2 (lower-left panel). In 
this case, both fit (EBL only and EBL plus ALPs) good. This is important, as it represents 
a good test of our ALP formalism: in those cases where the data are well understood and 
reproduced by conventional physics (regarding EBL and source modeling), ALP should not 
introduce any additional effect and should therefore be irrelevant if included. 

3.4 PG 1553+113 

This blazar was first detected in the VHE regime by the HESS collaboration in 2005 [48]. 
Since then, it has been extensively observed by the MAGIC collaboration, whose results are 
detailed in ref. [43]. In our study, the combined VHE data from this five-year campaign 
by MAGIC are used. The SSC modeling is also taken from ref. [43]. On the other hand, 
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the PG 1553+116 redshift has remained elusive (see e.g., ref. [49]), even though there are 
upper limits (z < 0.614 at la) coming from the expected EBL attenuation [49]. Recently, 
[50] estimated its redshift at z = 0.4 — 0.45 from spectral Ly-a signatures. For the redshift 
of the source we assume z = 0.4, warning the reader about its uncertainty. The result from 
applying our ALP formalism to this blazar is shown in figure 2 (lower-right panel). In this 
case, as in the PKS 1222+216 case, both the EBL-only and EBL-plus-ALPs curves provide a 
good fit to the observed data. Again, this tests that our ALP formalism does not introduce 
any additional effect in those cases where conventional physics by itself represents a good fit 
to the data. 

4 Discussion 

We have shown a possible alternative explanation to the pile-up problem observed at the 
highest energies in the EBL-corrected VHE spectra of the most distant blazars. This feature 
indeed sets a challenge to our current understanding of the emission mechanisms in blazars 
and/or the propagation of VHE photons through the IGM. While spectra corrected only for 
EBL absorption hardly describe some observations, the spectra corrected taking into account 
the existence of ALPs (with m a ~ 1 neV) provides a better fit to the data for all the objects 
discussed here, i.e. 3C 279, 3C 66A, PKS 1222+216 and PG 1553+113, independently of 
source modeling. Indeed, the pile-up problem is completely alleviated when including ALPs 
in the first two cases. On the other hand, in the case of PKS 1222+216 and PG 1553+113, 
the pile-up feature is not present and therefore ALPs are not needed; however ALPs provide 
a good fit as well, which actually represents a good check of our formalism (as ALPs would 
be always present if they exist). It is important to stress that these better fits to the data are 
not achieved by introducing any additional parameter in the fitting procedure. Indeed, the 
number of parameters involved in all the fits remains the same and equal to the degrees of 
freedom allowed by a simple power law. In this sense, we note that ALP parameters (such as 
E cr n) play a role which is just comparable to that played by those parameters describing the 
EBL model/physics. Both EBL-only and EBL+ALP curves shown in the plots correspond 
indeed to the outputs of different underlaying physical frameworks, each of one with their 
own parameters; both sets of parameters do not enter, in any case, in the posterior power-law 
fits to the data. 

It is also noticeable that those blazars that profit from the existence of ALPs roughly 
agree on the E cr n that better fits the data (~ 100 — 200 GeV). Indeed, for each particular 
case, the most appropriate E cr u lies near the lowest energy data point available. We note 
that every E cr a adopted was obtained by assuming two different Bjgm values; this fact 
points towards a degeneracy of our results for Bjgm- However, in practice, this degeneracy 
actually leads to very similar results (see Table 2). A way to break this degeneracy would be 
by observing AGN spectra to the highest possible energies, where different B values should 
lead to substantially different intensity boosts. 

We do not claim here to have shown that ALPs indeed exist. Other possible solutions to 
the pile-up problem within standard physics may exist such as large statistical and systematic 
uncertainties in the measurements, extensions of the SSC model [17, 19-22], or 7-ray photons 
from secondary interactions [17, 51]. Better source statistics, an extension of the analysis 
to lower-redshift sources, and a more detailed study are needed, but we aimed here to show 
explicitly in real VHE spectra of distant sources that the ALP possibility does deserve further 
study. We can envision multiple consequences of some of the ingredients with a prominent 
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role in the ALP formalism discussed here. For instance, note that in the case that ALPs 
exist with the properties adopted throughout this paper, the EBL upper limits derived from 
observation of VHE spectra of blazars [13-15] should be revised. 
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E mt [GeV] 


Source name 


no ALPs 


100 


200 


300 


400 


500 


3C 279 (2006) 


1.35 


0.37/0.42 


1.11/1.07 


1.10/1.07 


1.28/1.28 


1.36/1.36 


3C 279 (2007) 


3.75 


3.04/3.10 


1.96/2.19 


2.23/2.07 


3.74/3.74 


3.74/3.74 


3C 66A 


1.57 


0.18/0.24 


0.89/0.86 


0.88/0.87 


0.88/0.87 


1.34/1.34 


PKS 1222+216 


0.05 


1.47/1.35 


0.18/0.18 


0.09/0.14 


0.05/0.05 


0.05/0.05 


PG 1553+113 


0.43 


1.39/1.28 


0.19/0.18 


0.26/0.23 


0.32/0.32 


0.32/0.32 



Table 2. The \r values of the power-law fits to the de-absorbed spectra by the EBL effect only, and 
EBL plus ALPs for five benchmark values of the intergalactic E cr i t . Note that two Xr are given for 
every E crit , corresponding to the case with Biqm =0.1 nG and -B/gm = 1 nG, respectively. In bold, 
the best Xr f° r every source. Note that for PKS 1222+216 and PG 1553+113 cases, all E crit have 
roughly the same low Xri except E crit =100 GeV. 



- 11 - 



